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Abstract. We generalize the application of the functional renormalization group (fRG) method for the 
fermionic flow into the symmetry-broken phase to finite temperatures. We apply the scheme to the case 
of a broken discrete symmetry: the charge-density wave (CDW) mean-field model at half filling. We show 
how an arbitrarily small initial CDW order parameter starts to grow at the CDW instability and how 
it flows to the correct final value, suppressing the divergence of the effective interaction in the fRG flow. 
The effective interaction peaks at the instability and saturates at low energy scales or temperatures. The 
relation to the mean-field treatment, differences compared to the flow for a broken continuous symmetry, 
and the prospects of the new method are discussed. 

PACS. 71.45.Lr, 71.10.Fd Charge-density-wave systems, Lattice fermion models (Hubbard model, etc). 



1 Introduction 

Renormalization group (RG) techniques are a standard 
tool for the analysis of the low-energy physics of inter- 
acting fermions |1EJ2|. Although new effective degrees 
of freedom may be generated in the RG flow down to 
lower energies, perturbative RG flows with the initial de- 
grees of freedom - in our case electrons - are considered 
the least biased approaches to the variety of competing 
and conspiring tendencies at low scales. Without appro- 
priate self-energy corrections these flows often diverge at 
a nonzero energy scale, and not all fermionic modes can 
be integrated out. These divergences arc signatures of po- 
tential phase transitions, typically involving some kind of 
symmetry breaking. Hence, these flows still contain use- 
ful physical information. In many cases the physics below 
the critical scale can be explored by other means, e.g. by a 
mean-field treatment @| , by a bosonic description [5] or by 
exact diagonalization of a restricted low-energy Hamilto- 
nian [5]. Yet, the latter methods always require some sort 
of simplification of the low-energy action, and typically 
the final results depend on the scale where the transition 
from one treatment to the other is performed. Therefore, 
it is desirable to have an RG scheme which can be con- 
tinued beyond the noted divergence and which allows a 
controlled access to the low-energy phase. 

Recently, Salmhofer et al. have proposed an exten- 
sion of the fermionic functional RG flow into the symmetry- 
broken regime. The key idea is to include an arbitrarily 
small symmetry breaking field, which we include in the 
RG as a symmetry-breaking component in the initial con- 
dition for the fermionic self-energy. This component grows 



rapidly at the critical scale and prevents a true diver- 
gence of the interactions. Therefore, the flow can be con- 
tinued down to the lowest scales. Spontaneous symmetry 
breaking is captured by sending the external symmetry- 
breaking field to zero after integrating the flow. For mean- 
field models like the reduced BCS model, this procedure 
yields the exact solution already within a suitably chosen 
one-loop truncation of the vertex flow [7] . 

Of course this method should now be extended to other 
models and other physical situations. In particular, the 
hope is that non-mean-field models like the two-dimen- 
sional Hubbard model will be treatable. There, the method 
would yield the flow into symmetry-broken phases corre- 
sponding to correlations for which the initial interaction 
docs not contain an attractive component and where the 
attraction is generated during the flow. For example, this 
is the case for <i-wave pairing in the Hubbard model on the 
two-dimensional square lattice near half filling. Further- 
more, the non-mean-field parts of the interactions would 
not need to be dropped and their influence could be stud- 
ied. 

Before proceeding to such more complicated problems, 
we extend the application of the formalism to nonzero 
temperatures in the context of the charge-density wave 
(CDW) mean-field model on a d-dimensional hypercubic 
lattice at half band filling. Again, we closely monitor the 
correspondence to mean-field theory. This example is in- 
teresting because, in contrast to the BCS model studied 
in Ref. only a discrete symmetry is broken. Hence we 
learn about the differences between the RG flows for bro- 
ken continuous and discrete symmetries. 



2 



Gersch, Honerkamp, Rohe, Metzner: RG flow into CDW state 



This paper is organized as follows. In Sec. [3 we intro- 
duce the CDW model which we analyze in the remainder. 
In Sec. |3| we present an exact solution in the thermody- 
namic limit employing a rcsummation of perturbation the- 
ory. We briefly introduce the renormalization group tech- 
nique in Sec. 0] and write down the functional RG equa- 
tions for the self-energy and the effective interaction. In 
Sec. [3 these equations are discussed in the case of vanish- 
ing temperature, where they can be presented in a consid- 
erably simpler way. The finite-temperature equations are 
treated numerically in Sec.[S] Conclusions and outlook are 
presented in Sec. 



2 Model 

We consider spinless fermions on a cf-dimensional hyper- 
cubic lattice with TV sites labeled by x. The kinetic energy 
shall be given by a tight-binding dispersion with nearest- 
neighbor hopping amplitude t. Further we assume a repul- 
sive nearest- neighbor density-density interaction Vq. Then, 
the Hamiltonian reads 

x,n x,n 

The sum over n runs over the d unit vectors of the d- 
dimcnsional hypercubic lattice with lattice constant set 
to one. Now we assume half filling with an average of one 
particle per two lattice sites, (n x ) = 1/2. In this case it 
is easy to see that the nearest-neighbor density-density 
repulsion will favor a charge density wave (CDW), i.e. a 
periodic arrangement of the particles in which the proba- 
bility of two particles being nearest neighbors is reduced. 
This tendency competes with the hopping term of the ki- 
netic energy and with the entropy which arc minimized 
and maximized, respectively, for an equal population of all 
sites. In the case where a CDW is formed with a fixed mod- 
ulation amplitude ncDW, the charge density (n x ) takes 
only two values, l/2±ncDW; on the two inter-penetrating 
sublattices with doubled unit cells. Therefore, the half- 
filled CDW state breaks a discrete Ising-type symmetry, 
and the degeneracy of the ordered state is twofold. While 
in one dimension such order can only occur in the ground 
state, in two and higher dimension it is possible to have a 
nonzero critical temperature T c . 

For general band filling, a CDW with an incommen- 
surate modulation wave-vector generates infinitely many 
different values of average populations on the lattice sites. 
This corresponds to an arbitrary phase of the density mod- 
ulation at an arbitrary reference lattice site. Hence in the 
incommensurate case the CDW breaks a continuous sym- 
metry. 

Continuing with the half-filled case, we Fourier-trans- 
form the Hamiltonian using 

^ = ^E eife % ( 2 ) 



and obtain 

ff = -2t£i(fc)cta 

k 

+ 77 £ (9) c * 1 c fc 1 -, c L c * a +,- ( 3 ) 

k± ,k2 ,q 

Here we have introduced L(k) = Yli=i cos(fcj). Abbreviat- 
ing the dispersion as £(fe) = — 2t Yli=i cos(fci) = —2tL(k), 
we note that the nesting condition £(fc) = — £(fc + Q) is 
fulfilled for the wave- vector Q = (tt, . . . ,vr). This causes a 
divergence of the non-interacting charge response at the 
nesting wave- vector for T — ► 0. If the interactions are 
treated in the random phase approximation, the charge 
response at Q will actually diverge at a finite temperature, 
giving an estimate of T c for the charge-density wave forma- 
tion. In a mean-field treatment of the symmetry-breaking 
in which the interaction term is decoupled with an al- 
ternating charge density, anomalous particle-hole pairing 
expectation values {c\ + qC^} will become nonzero below 
the transition temperature. 

Although this physical picture of a CDW transition is 
essentially correct regarding the ground state properties, 
the model (j2J) cannot be solved exactly. This changes if 
we reduce the interactions by only keeping processes that 
change both particle's momenta by Q, 

fc 

E C k 1 C k 1 -Q C k 2 C k 2 +Q- ( 4 ) 
fci ,k 2 

Now the interaction between fermions on the lattice sites 
x and x' corresponds to an infinite range density-density 
interaction with oscillating sign, 
— VoiV -1 cos[Q(a; — x')]n x n x i. 

Analogous to the reduced BCS pairing model [S], in 
the reduced model J3J mean-field theory becomes exact 
in the thermodynamic limit N — > oo and there is a CDW 
ordered state below a critical T c > in any dimension 
d. For the half-filled commensurate case, the CDW state 
is twofold degenerate and the electronic order parameter 
ncDW is real. 

Next we include a real external field A ext cos(Qx). 
This term breaks the translational symmetry of the orig- 
inal Hamiltonian explicitly and lifts the degeneracy of 
the two CDW ground states. In momentum space this 
term couples to pairs of particles whose momenta differ by 
Q = (7r, . . . , 7r). If Z\ ox t is very small compared to all other 
relevant energy scales, it does not change macroscopic ob- 
servables away from the critical temperature. However, it 
allows us to integrate the RG differential equations over 
all scales without encountering divergences. The Hamilto- 
nian including the external field reads 

H rc d = -Skin + ^2 ^ cxtC fc c fc+Q 

k 

_ JJ C k l c k l +Q c k 2 C k 2 -Q- (^) 

fci ,k 2 
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Wc introduce a frequency-space field-integral represen- 
tation in the usual way. Writing T for temperature, we de- 
fine the fermionic Matsubara frequencies uij := (2j + 1)ttT 
and obtain Grassmann fields ipk,u n > ^'fc.^n • Introducing the 
Nambu-likc notation 



the partition function reads 



D(tj},ijj) exp 



n.h 



T 3 



fei ,k 2 



(6) 

Here, lu^ = lj„ 1+ „ 2 _„ 3 and Q is the matrix inverse of the 
free Green's function (see e.g. [§]) 



G (£,w n ) 



1 



-lU) n 

-A, 



(7) 



We have dropped the momentum argument of £ for brevity. 
The minus sign in front of the first £ arises from the Q- 
anti-pcriodicity of the cosine. The doubling of the number 
of degrees of freedom by introducing the Nambu nota- 
tion is compensated by the 1/2 in front of the quadratic 
part of the action. Alternatively, we could integrate over 
only half the Brillouin zone. Due to the explicit symmetry 
breaking, there is a nonzero off-diagonal component in the 
propagator. If we apply mean-field theory by decoupling 
the interaction term in the CDW channel, the CDW order 
parameter due to the interaction, A m , will add to Zi e xt- 



3 Mean-Field Solution and Resummation 

In this section, we resum the perturbation expansions for 
the self-energy and the effective interaction in the thermo- 
dynamic limit N — > oo. For our model this is equivalent to 
mean-field theory in the CDW amplitude with modulation 
wavevector Q. 

Due to the explicit symmetry breaking with Z\ cxt 7^ 
there is a nonzero frequency independent off-diagonal self- 
energy which scatters a fcrmion with wavevector k to 
k + Q and vice versa. Therefore, besides normal (diago- 
nal) wave-vector conserving propagators, also anomalous 
(off-diagonal) propagators appear in the diagrams of the 
perturbation theory. As indicated above, compared to the 
non- interacting propagator in JJJ), the effect of the inter- 
actions is that we have to supplement Zi e xt with an off- 
diagonal self-energy Z\; a caused by the interaction. In the 
following we determine its value. 

Since the interaction term in JSJ includes a prefactor 
1 /N, only diagrams having one summation over all degrees 



(a) 




(b) 



Fig. 1. Resummation of the perturbation expansion of the 
reduced charge-density-wave model: (a) self-energy (hatched 
circle), and (b) effective interaction (hatched rectangle). Undi- 
rected internal lines carry a summation over their Nambu-like 
indices in addition to the frequency summation and momentum 
integration. Double-arrowed lines correspond to off-diagonal 
propagator elements, bold lines to full and thin lines to bare 
propagators. 



of freedom per interaction line contribute in the thermo- 
dynamic limit. This together with the special structure 
of the interaction entails that only the bubble diagrams 
shown in Fig. 1(b) contribute to the perturbation series 
for the effective interaction in second order of the bare in- 
teraction. In all other second order diagrams the restricted 
interaction implies that the internal momentum is fixed by 
an external momentum. Then the contribution of these di- 
agrams is 0(1/N) and vanishes for N — > 00. Furthermore, 
since each bubble in the remaining diagrams transfers a 
momentum of Q due to the special structure of the bare 
interaction, only effective interaction processes with a mo- 
mentum transfer of Q exist. This causes the normal diag- 
onal self-energy to vanish as well. These arguments can 
be iterated to arbitrary order, and hence only the bubble 
chains indicated in Fig. 1(b) contribute to the perturba- 
tion expansion of the effective interaction. In particular, 
no anomalous effective interactions which violate the orig- 
inal translation symmetry are generated. This is different 
away from half-filling where pairs of internal anomalous 
propagators can generate effective interactions which vi- 
olate momentum conservation by adding even multiples 
of the modulation wavevector Q to their incoming total 
momentum. In the half-filled case 2Q = (2tt, . . . , 27r) is 
a reciprocal lattice vector, and therefore no new effective 
interactions appear. 

In Fig.^a) we show the perturbation expansion for the 
off-diagonal self-energy. Again, due to the special structure 
of the interaction only those diagrams which retain one 
momentum integration per interaction line contribute to 
the perturbation expansion of the self-energy in the ther- 
modynamic limit. Furthermore, the interaction is nonzero 
only for momentum transfer Q. Thus, in the thermody- 
namic limit only diagrams constructed from tadpole di- 
agrams contribute. If a tadpole diagram is connected to 
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at most one other tadpole diagram, only the terms off- 
diagonal in the Nambu-like space are non-vanishing. This 
is because the outgoing momentum of such an off-diagonal 
line equals the incoming momentum plus Q. The exact 
resummation is equivalent to self-consistent Hartree-Fock 
theory of which the Fock contributions vanish in the ther- 
modynamic limit. Another way of saying this is that Har- 
tree mean-field theory is exact in this model. 

In the following, we assume the density of states to 
be constant over the Brillouin zone and equal to its value 
at the Fermi energy to remove dimension-specific effects. 
Evaluating Fig. Ufa) and denoting the self-energy due to 
the interaction by A- m: the band edge by W, the density 
of states at the Fermi energy po, £ :— A- la + A cxt as well 
as E := \/£ 2 + £ 2 , we obtain the gap equation (compare 
e.g. PU) 



0.2 



A r 



V £p 



w tanh(£/2T) 
^ E 



OS) 



which resembles the BCS ^T] and excitonic insulator ^21 
gap equations. If we set T = A ext = 0, solving JSJ analyt- 
ically yields (compare e.g. (11113) ) 



2t- 



1 



sinb.(l/VoA)) 



V po«l , / 1 

4icxp - 



VqPq 



Similar equations are found in the BCS theory for super- 
conductors [TT] . 

Wc now turn to the effective interaction with zero fre- 
quency transfer. As argued above, only fcrmionic bubbles 
contribute and all other diagrams vanish in the thermo- 
dynamic limit N — > oo. In the notation, we omit the ar- 
guments of the hyperbolic functions originating from the 
Fermi distribution from now on: they are always E/2T. 
Abbreviating the bubble integral 



; w d£ 



£ 2 -2 r 

— cosh H tanh 

2T E 



the evaluation of Fig. QJb) yields 

V 



V 



1-V B 



(10) 



(11) 



If A cxt = 0, we can define the critical temperature T c via 
the condition that 1 = VoB(T c ), i.e. that the denominator 
of 1)11( 1 vanishes and V diverges. By analyzing 1 )10(1 we note 
that this criterion for T c coincides with the temperature 
below which we can find a nonzero solution of the gap 
equation (JSJ). For small Vq and £ = wc find 1.76T C ~ 
4iexp(-l/Vb/Oo) ~ S(T = 0) (see e.g. [El). This is true 
when the approximation in ||5J) is valid. 

As a finite Z\ ex t will later allow us to integrate the 
renormalization group equations over all scales, we are in- 
terested in the dependence of the mean-field solutions JSJ 
and 1)11(1 on Z\ oxt . In particular, we will illustrate where 
the system including a finite symmetry-breaking field 
is practically indistinguishable from the system without 
symmetry-breaking field, To this end, we show in the 
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Fig. 2. (top) Solutions S(T) of the gap equation for Vo = 2t 
and small to intermediate initial self-energies zA ex t- (bottom) 
Temperature dependence of the effective interactions calcu- 
lated by resumming the perturbation expansion for Vo — It 
and small to intermediate initial self-energy, po = l/2n for all 
numerical calculations. 



upper part of Fig. convergence of £ A <=^ (T) to £(T) 
for Z\ cx t — ► 0. The convergence is worst in the vicin- 
ity of a kink that can be discerned at T w 0. It in the 
graph for Z\ ext = 0. We would like to be more concrete 
about the temperature range in which the system includ- 
ing a symmetry-breaking field can be used to approximate 
the system without symmetry-breaking field. To deter- 
mine more precisely the temperature range in which the 
symmetry-breaking field has little impact on the physics 
of the system, we consider the dependence of the effec- 
tive interaction on both temperature and the symmetry- 
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breaking field as plotted in the lower part of Fig. Its 
singularity is regularized by a nonzero A cxt which cuts 
off the integrand of the bubble integral 1|10|) . The figure 
also shows a strong suppression of the effective interaction 
peak with increasing Z\ oxt (the dependence of the peak 
height on A ext is plotted in Fig. [TJ) . Its location converges 
to T c from above for Z\ cx t — + 0. In the figure, we also spot 
convergence of V A "^ (T) to V(T) for T < T c and T > T c . 
We see that V A ^(T) « V(T) if T is outside the double 
width at half maximum of V A ~*(T). If T is outside said 
width and also below T c , we learn from the upper part of 
the figure that E A <** (T) w E(T). This is thus the temper- 
ature region in which the initial symmetry-breaking field 
does not appreciably change the physics. 



4 Renormalization Group Setup 

We will now apply the one-particle irreducible (IPI) func- 
tional renormalization group scheme to the charge-density- 
wave. The scheme is described in pQ (see also Q3]) and is 
here employed in the version suggested by Katanin ^H] • In 
this version, the differential equations for the self-energy 
and four-point function constitute a closed system. The 
special structure of our bare interaction will allow us to 
further simplify the equations. We will sec that in contrast 
to the BCS flow no anomalous effective interactions 
are generated. In later sections, we will verify numerically 
that the fRG reproduces the resummation results. 

Following pQ, we introduce a positive real parameter 
s and a cutoff A := 2iexp(— s). In our calculations, all 
modes satisfying |£| > A have been integrated out. The 
effective interaction and self-energy at this scale can thus 
be interpreted as parameters of an effective theory for 
a reduced system with smaller bandwidth. We start the 
flow at s = (A = 2t), where all modes have yet to be 
integrated out, and integrate the flow down to A = 0. 
To analytically implement this procedure, we introduce 
x(£, A) as a placeholder for any cutoff function, which shall 
have range [0, 1], assume the value 1/2 when A = £, ap- 
proach when reducing £ below A and 1 when increasing 
£ above A. We replace Q(£) by Q(£)/x(£, A), suppress- 
ing low-energy modes in the field integral and rendering 
the self-energy and effective interaction scale-dependent. 
In the following, dots over symbols will denote differenti- 
ation with respect to s. In the diagrammatic representa- 
tions, a heavy line represents a scale-dependent full prop- 
agator (calculated with the scale-dependent self-energy), 
a hatched circle represents the scale-dependent self-energy 
and a hatched rectangle represents the scale-dependent ef- 
fective interaction. We furthermore introduce the single- 
scale propagator S := GQG (all scale-dependent). A line 
drawn through propagator lines indicates an s-derivative. 

In the flow of the effective interaction, having d s (GG) 
[T5] instead of §G + G§ as in the original IPI RG scheme 
PP is crucial in order to be able to follow the flow down 
to A = 0. In the Katanin version of the IPI scheme, the 
exact hierarchy of RG differential equations with the ex- 
ception of the first equation is written using only full four- 
point functions, full propagators and scale-differentiated 



full propagators on the right-hand sides. The differential 
equations for the self-energy and the four-point function 
thus constitute a closed system. In particular, the con- 
tribution to the flow of the four-point function which in- 
volves the six-point function in the standard IPI scheme 
PP is taken into account in the Katanin scheme by dia- 
grams involving only four-point functions and full or scale- 
differentiated full propagators, respectively. The part with 
overlapping loops of this contribution vanishes for the spe- 
cial structure of the interaction and the thermodynamic 
limit considered here. The remaining part of this contribu- 
tion is taken into account by the replacement of the single- 
scale propagator by the scale-differentiated full propaga- 
tor. Hence, the graphical form of the fRG differential equa- 
tions is as shown in Fig. 

In the flow of the effective interaction (right part of 
Fig. m, the summation over the 'Nambu' indices of the 
internal lines includes normal and anomalous propaga- 
tors. Analogous to the perturbation expansion in the pre- 
vious section, in the limit N — > oo the special momentum 
structure of the initial interaction is conserved by the RG 
flow of the effective interaction. In the present half-filled 
case, no new effective interactions with different external 
legs are generated. This is in contrast to the BCS pairing 
model, where U(l)-symmetry-breaking vertices with four 
incoming or four outgoing legs are generated. Linear com- 
binations of normal and anomalous components of the ef- 
fective interacion could be identified with 'amplitude' and 
'phase' vertices. In the half-filled CDW model the effective 
interaction behaves like the 'amplitude' vertex of the BCS 
problem, while the 'phase' vertex is missing as there is no 
breaking of a continuous symmetry. The self-energy only 
contains the off-diagonal part of the single-scale propaga- 
tor, again due to the restricted form of the interaction. 
For a derivation see pQ and jZj. Ref. also provides an 
analytical proof of the equivalence of the flow equations 
(Fig. |2J and the resummation equations (Fig. [JJ. 



Fig. 3. The hierarchy of flow equations for the self-energy and 
the effective interaction. 



Evaluating the diagrams of Fig. [31 we suppress scale 
and frequency dependence in the notation and write E := 
V £ 2 + X 2 ^ 2 . The RG differential equations now read 



A ( A 2 
V= PQ V Z {- -^tanh 



3^ 
E 2 



E \E 2 



d£ EE 



— cosh 2 

2E 



tanh — — — cosh 2 
2 



-P 2 E< 



tanh 



2 cosh^ 



(12) 
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E 



(13) 



where we have used § = xdG/dx, the limit of a sharp 
cutoff, and (3.19) of (The simplification of the equa- 
tions obtained using the limit of a sharp cutoff is due to 
Tilman Enss). The initial conditions arc ^l^o = A cxt 
and V | s= o = Vo. If we choose Z\ e xt = 0, £ remains zero 
for all scales. Then only the first term in Q12f) remains, and 
the flow will diverge at a nonzero scale A c . We would thus 
be neither able to integrate over all scales nor compare 
our results to the resummation of section [3] We therefore 
always study the RG equations for a finite A CKt , usually 
10~H. We will show that we can thereby circumvent the 
divergence of the effective interaction and still get arbi- 
trarily close to the exact mean-field results without ex- 
plicit symmetry breaking in the temperature range deter- 
mined in section [3] 



5 Renormalization Group at Zero 
Temperature 

We first consider zero temperature, T = 0. In this case, 
it is easier to write down the analytical expressions, the 
roles of the different right-hand side terms are easily iden- 
tifiable, and the value of the effective interaction at A = 
is nearly independent of Z\ ox t, as A c is larger than the 
double width at half maximum of the effective interaction 
flow peak (see Fig. 0] lower part). 

In the limit T -> 0, (JT2J) and JT3J) become 



.±y=y* A2po 
dA VA^~ 



3V 2 £ 



£2- 



d£ 
dA 



d£ 



~£ =V£p -==L . 

dA m v^ 2 + £ 2 



epo 

(14) 
(15) 



Notice that we have replaced the usual derivative with 
respect to s by —Ad/dA. The typical shapes of this flow 
are exhibited in Fig. We see that the divergence of the 
effective interaction is regularized by the initial symmetry- 
breaking field. Hence, we can integrate out all modes of the 
fcrmionic spectrum. The four-point function in the CDW 
problem behaves analogously to the linear combination 
of normal and anomalous four-point functions in the BCS 
problem 7 which drives the flow of the pairing amplitude. 
The yl-dependence of £ around the scale where V peaks 
approaches a kink for zi cxt — > 0. The graph of £(A) re- 
sembles the graph of £{T), showing that temperature acts 
in a similar fashion as the cutoff in this system. The main 
difference between the graphs is that £(T — > 0) saturates 
exponentially while £{A — > 0) is linear. 

The impact of the self-energy feedback on the RG flow 
can be thoroughly understood by analyzing the terms ap- 
pearing in (|14|) and l|15(l . Bear in mind that we think of the 
flow as progressing from larger to smaller values of A, i.e. 
from right to left in our plots. We first consider the limit 
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Fig. 4. Low-energy portions of the T — flows of the self- 
energy (top) and effective interaction (bottom) for Vo — 2t 
and various zi cx t. po — l/2n for all numerical calculations. 
Compare with Fig. [5] 



£ = 0. Now, the self-energy no longer flows while the effec- 
tive interaction flows according to —dV/dA oc V 2 / A. For 
any positive Vo, the solution of the corresponding initial 
value problem is singular, showing that at T = symme- 
try is broken for all repulsive initial interactions. 

An arbitrarily small initial symmetry-breaking field 
immediately has dramatic consequences: As soon as V 
begins to increase strongly, so does the self-energy due 
to the coupled effects of a large V and a back-feeding £ 
on the RHS of |(T5J|. The second term on the RHS of (fUjl 
is negative. It corresponds to a correction of the flow by 
the modes which have been integrated over at A' = £ > A 
reflecting that their self-energy £ is changing in the flow. 
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It becomes large as the slope of £ exceeds £ itself, while 
the positive first term is damped by an effective A 2 /£ 3 - 
dependence as soon as £ becomes much larger than A. 
The effective interaction is hence pulled back down. It 
never reaches zero since the negative second term on the 
RHS of l|14|) becomes proportional to V 3 while the pos- 
itive first term becomes proportional to V 2 (d£/dA is 
0(V)). This implies that the self-energy never decreases. 
The two terms on the RHS of (|f 4fl thus approach mutual 
cancellation, causing V and in turn —d£/dA to be almost 
constant in A (the contributions from £ on the RHS of 
(JEJ) cancel when £ > A). 

The behavior described above remains qualitatively 
the same at finite temperatures, as long as these are below 
the double width at half maximum of V A ° xt (T) . This is 
illustrated by the lowest-temperature curve of the lower 
panel of Fig. El The shape of V(A) is very similar to the 
shape of V{T) in the lower part of Fig. El As it is found 
in many one-dimensional problems, temperature thus has 
an effect comparable to that of an energy or momentum 
cutoff, respectively, as implemented here. 

6 Renormalization Group at Finite 
Temperatures 

We now turn to the analysis of the finite temperature fRG 
equations (|I2J) and (|I3Jl . Due to their more involved na- 
ture, the analysis is largely numeric. We find the same 
behavior as calculated using the conventional resumma- 
tion methods applied in section El The initial symmetry- 
breaking field plays the same crucial role as in the zero- 
temperature case. We use a coupling of Vq = 2t in the 
following. 

The upper plot of Fig. El shows how the flow of the 
self energy flattens out when the temperature is increased 
beyond T ks O.li. Above T c , £{A = 0) vanishes in the 
limit A cxt — > 0. The RG flow of the effective interaction 
is shown in the lower panel of Fig. El The graph of the 
flow is pushed to the left with increasing temperature, 
its shape remaining largely unchanged. As the maximum 
approaches zero, the final value of the effective interaction Fig. 5. Flows of the self-energy E and the effective interac- 
increases until the maximum has reached A = 0. This tion V plotted against A for temperatures around the transi- 
corresponds to the behavior of the effective interaction tion, Vq = 2t and A ext = 10~ 4 t. p — l/2n for all numerical 
on the low-temperature side of Fig. El (lower part). For calculations, 
even higher temperatures, the final value of the effective 
interaction decreases, corresponding to the behavior on 

the high-energy side of Fig. El zer0 sc& l e - However, the upper part of Fig. gl implies that 

The lower part of Fig. El shows that A C (T) saturates lim^^o £a^{A) exists and is approached in a contin- 
quickly far below the critical temperature. It also illus- uous fashion. 

trates the motion of the effective interaction flow maxi- For T > 0, if zA ext is one or two percent of the final £, 
mum, which is pushed to lower scales by an increase in there is a clearly observable steep rise (see the upper part 
temperature as described above. It approaches zero in a of Fig. EJ which would allow a rather precise determination 
linear fashion with increasing temperature, in contrast to of T c even in more complicated models. However, note that 
£(T — > T c ), which exhibits a square-root behavior. for these values of A eyLt , the effective interaction at the 

We are furthermore interested in the dependence of critical scale still reaches values ~ lOOi which are much 
the flow of the self-energy at finite temperatures on Zi e xt- larger than the bandwidth At (see, e.g., the lower part 
In contrast to the resummation treatment case (see the of Fig. El for the temperature evolution) . The agreement 
upper part of Fig. EJ, we cannot set Z\ e xt to as the ef- between fRG and mean-field results despite the fact that 
fective interaction would diverge before the flow reaches the effective interaction has grown to values large enough 
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Fig. 6. (top) Self-energy against temperature calculated us- 
ing fRG methods for Vo = 2t and small to intermediate A e *t- 
Minute deviations in comparison to Fig.|5]are due to numerical 
errors in the integration of the fRG differential equations, (bot- 
tom) Position of the maximum of the cutoff-dependent effective 
interaction against temperature for small to intermediate Ziext- 
po = 1/2-7T for all numerical calculations. 



to render any perturbative scheme useless in any general 
model underlines that the truncated fRG is exact for our 
model. 

At the T-dependent critical scale A c , the effective in- 
teraction reaches a maximum whose height depends sin- 
gularly on Z\ ex t (see Fig. 0). A numerical analysis shows 
that the maximal effective interaction is oc 1/Z\" xt with 
a ~ 2/3. The values for the maximal effective interaction 
can be read off for T = in Fig. 0] (bottom) and for finite 
temperatures in the lower part of Fig.0 For a straightfor- 



ward application of the method to non-mean-field models, 
it may be necessary to restrict the maximum of the effec- 
tive interaction to smaller values to justify a perturbative 
treatment. Then the rise of Sa=o(T) at T c gets smeared 
out strongly (see Fig. (upper part) and Fig. (upper 
part)). The upper part of Fig.j3sh.ows that far below T Ci 
the relative error AS / E is linear in Zi ox t. At the high- 
est considered Z\ ex t, it is already of the order of 0.1 while 
V(A C ) is still larger than two times the bandwidth. 
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Fig. 7. (top) Relative deviation of the fRG-derived self-energy 
from the value without external field against the initial anoma- 
lous self-energy for temperatures from T = to just above T c . 
The values for T = 0.12f are 1 and thus not visible in the 
selected range, (bottom) Maximal value of the effective inter- 
action during the flow against the initial anomalous self-energy 
for temperatures from zero to just above T c . 
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7 Conclusion 

Wc have applied the functional renormalization group to a 
simple model for the CDW transition of spinlcss fermions 
on a hypercubic lattice in d spatial dimensions. The main 
purpose was to test and further understand the recently 
developed method for the fRG flow into symmetry-broken 
phases in the case of a broken discrete symmetry and to 
generalize it to finite temperatures. 

Earlier RG treatments of low-dimensional interacting 
fermions generically encountered flows to strong coupling 
where a class of coupling constants becomes too large to be 
treatable pcrturbatively. Then, the flow had to be stopped. 
While these flows still contain ample information about 
the low-energy physics of the system, a continuation down 
to lowest scales is desirable to obtain a controlled access 
to possible symmetry-broken phases. A recent treatment 
of the BCS model showed that the divergence of the inter- 
actions at finite RG scale can be avoided by introducing 
a small initial symmetry-breaking field. Furthermore, a 
correct treatment of the flow of the corresponding anoma- 
lous self-energy yields the correct value for the energy gap 
at zero scale. The interaction vertex corresponding to the 
Goldstone boson in the case of breaking of the continuous 
U(l) symmetry still becomes large oc 1/Z\ cxt when the 
initial symmetry-breaking field Z\ e xt is sent to zero. 

Here we applied this method to the CDW problem 
with a half-filled band. Including a small initial anoma- 
lous CDW self-energy A oxt / we found that the flow of 
the effective interaction goes through a maximum at the 
critical scale A c and falls off below. In this case, the effec- 
tive interaction does not become large as A — > and no 
anomalous interactions are generated, reflecting the fact 
that the discrete symmetry breaking of the commensurate 
CDW does not produce a Goldstone mode. At A c , the 
anomalous self-energy U corresponding to the strength of 
the discrete symmetry breaking rises steeply from its ini- 
tial value A cxt and finally approaches a finite value for 
A — > 0. Letting Z\ ox t — > 0, we found that Ea=o converges 
to the exact result from the case without explicit symme- 
try breaking, Z\ ext = 0. 

Studying the flow at nonzero temperatures, we have 
noted that the temperature dependence of the critical 
scale differs qualitatively from the temperature depen- 
dence of the off diagonal self-energy, i.e. no square-root 
like rise of A C (T) is found below T c . Therefore, A C (T) is 
only a rough estimate for the self-energy. Thus, a proper 
treatment of the flow of the self-energy into the symmetry- 
broken phase is needed in order to assess the temperature 
dependence of the spectral gap. 

Similar to the BCS problem, the maximal value of the 
effective interaction at A c depends crucially on the ex- 
plicit symmetry breaking Z\ cx t. For the mean-field model 
we have investigated in the present work, the truncated 
RG scheme is exact. Thus, large effective interactions do 
not cause problems. However, for the practical applica- 
tion of the fRG scheme to non-mcan-field models where 
our truncated RG scheme is not exact it is necessary that 
the interactions do not become too large at any scale. 
Otherwise the feedback of the interactions in the ordering 



channel on processes in other channels would become large 
and the perturbative RG would break down. In the case 
of breaking of a discrete symmetry studied here, the only 
region where the interactions become large at zero tem- 
perature is near A ~ A c . There, however, our explicit cal- 
culation shows that if we want to stay close to the mean- 
field results without explicit symmetry breaking wc cannot 
choose too large a Z\ ox t- This in turn leads to maximal val- 
ues of the effective interaction which are much larger than 
the bandwidth. If one wants to work with couplings which 
stay within the order of the bandwidth (or expressed dif- 
ferently, with dimensionless coupling less than one), the 
deviations from the A ext = results are rather large. At fi- 
nite temperatures, the explicit symmetry breaking smears 
out the transition. Again, the onset of the growth of the 
off-diagonal self-energy only allows for a reasonable defini- 
tion of a critical temperature when the maximal effective 
interaction is allowed to get large. 

On the one hand, the large effective interactions dis- 
cussed above could spoil the application of this truncated 
RG scheme to non-mean-field models. On the other hand, 
one should keep in mind that in non-mean-field models 
large effective interactions only occur for a narrow range of 
frequencies and momenta, if the bare interaction is weak. 
Hence, the feedback of these large interactions into other 
interaction channels will be strongly reduced by phase- 
space factors. It will be interesting to see how this physics 
saves or destroys for more general types of interactions 
the mean-field picture of the reduced models. In particu- 
lar, the space dimension plays a crucial role in this respect. 
Understanding these issues could provide insights into the 
emergence of short-range ordered states in a weak cou- 
pling picture. 

With the steps outlined in this article, we believe that 
the fRG flow into symmetry-broken states is well-under- 
stood for mean-field- type models. In the future it can be 
generalized to more complicated models, where it will re- 
veal its full usefulness. Finally, wc note that an alternative 
method for the flow into the symmetry-broken phase based 
on a two-particle irreducible RG formalism has been put 
forward by Dupuis |17| . It will be interesting to explore 
the applicability of these new RG schemes in more realistic 
models. 

We thank Tilman Enss, Andrey Katanin, Julius Reiss 
and Manfred Salmhofcr for useful discussions. 
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